The Spitzer Survey of Interstellar Clouds in the Gould Belt. V. 
Ophiuchus North Observed with IRAC and MIPS. 



J. Hatchell 

Astrophysics Group, Physics, University of Exeter, Exeter, EX4 4QL, United Kingdom 

hat chell@astro . ex . ac . uk 

S. Terebey 

Department of Physics and Astronomy PS315, 5151 State University Drive, California State University at 

Los Angeles, Los Angeles, CA 90032, USA 

sterebeOcalstatela. edu 

T. Huard 

Department of Astronomy, University of Maryland, College Park, MD 20742, USA; Columbia Astrophysics 
Laboratory, Columbia University, New York, NY 10027, USA 

thuardOastro . umd . edu 

E. E. Mamajek 

University of Rochester, Department of Physics & Astronomy, P.O. Box 270171, Rochester, NY, 

14627-0171 USA 1 

emamaj ekOpas . rochester . edu 

L. Allen 

National Optical Astronomy Observatory, 950 N. Cherry Ave., Tucson, AZ 85719, USA 

lallenOnoao . edu 

T.L. Bourke 

Harvard- Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA 
tbourke@cf a . harvard . edu 

M. M. Dunham 

Department of Astronomy, Yale University, Box 208101, New Haven, CT 06520-8101, USA 

michael . dunhamOyale . edu 

R. Gutermuth 

Department of Astronomy, University of Massachusetts, Amherst, MA 01002, USA 

rgutermu@astro . umass . edu 



1 



P. M. Harvey 

Astronomy Dept., University of Texas, Austin, TX 78712, USA 
pmMastro . as . utexas . edu 



J. K. J0rgensen 

Niels Bohr Institute and Centre for Star and Planet Formation, University of Copenhagen, Juliane Maries 

Vej 30, DK-2100 Copenhagen 0, Denmark. 

jeskj Onbi.dk 

B. Merin 

Herschel Science Centre, European Space Astronomy Centre (ESA), P.O. Box, 78, 28691 Villanueva de la 

Canada, Madrid, Spain 

Bruno .Merin@scd.ops . esa. int 

A. Noriega-Crespo 

Science Center, MC 220-6, California Institute of Technology, Pasadena, CA 91125, USA 

albertoOipac . caltech . edu 
and 



D. E. Peterson 

Harvard- Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA 
dpetersonOcf a. harvard. edu 



ABSTRACT 

We present Spitzer IRAC (2.1 sq. deg.) and MIPS (6.5 sq. deg.) observations of star formation 
in the Ophiuchus North molecular clouds. This fragmentary cloud complex lies on the edge of 
the Sco-Cen OB association, several degrees to the north of the well-known p Oph star-forming 
region, at an approximate distance of 130 pc. The Ophiuchus North clouds were mapped as part 
of the Spitzer Gould Belt project under the working name 'Scorpius'. In the regions mapped, 
selected to encompass all the cloud with visual extinction Ay > 3, eleven Young Stellar Object 
(YSO) candidates are identified, eight from IRAC/MIPS colour-based selecti on and three from 



2MASS Ks/MIPS colours. Adding to one source previously identified in L43 (jChen et al.ll2009f) . 
this increases the number of YSOcs identified in Oph N to twelve. During the selection process, 
four colour-based YSO candidates were rejected as probable AGB stars and one as a known 
galaxy. The sources span the full range of YSOc classifications from Class 0/1 to Class III, and 
starless cores are also present. Twelve high-extinction (Ay > 10) cores are identified with a total 
mass of ~ 100 Mq. These results confirm that there is little ongoing star formation in this region 
(instantaneous star formation efficiency < 0.34%) and that the bottleneck lies in the formation 
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of dense cores. The influence of the nearby Upper Sco OB association, including the 09V star 
C Oph, is seen in dynamical interactions and raised dust temperatures but has not enhanced 
levels of star formation in Ophiuchus North. 

Subject headings: circumstellar matter - infrared: stars - stars: formation - stars: evolution 



1. Introduction 

The Ophiuchus North (Oph N) molecular 
clouds lie 20° above the Galactic plane in the 
direction of the Galactic Centre. They are part of 
the same filamentary cloud complex as the well- 
studied Ophiuchus L1688 and L1689 clouds, but 
lie several degrees to the north, on the boundary 
of the constellation of Ophiuchus with Scorpius. 
Like the Ophiuchus L1688 and L1689 clouds, 
they are illuminated from the northwest by the 
Upper Scorpius subgroup of the Sco-Cen OB as- 
sociation. 

The region has been st udied little. Early 
low-reso l ution CO mapping (|de Geus et al.lll990t 
de Geusl Il992l) showed the filamentary structure 
of the clouds in the Ophiuchus regi on (mirrored in 



the e xtinc tion maps published bv | Dobashi et al 



20051 and iRowles fc Froebrichl 120091 ) . and sug- 
gested a shock origin due to expanding shells sur- 
rounding the Upper Sco subgroup. A detailed 
stu dy of the molecular clouds was made in 13 CO 
by iNozawa et al. ( 1991 ) which gives an excellent 
overview of the cloud complex, and its relation- 
ship to the Ophiuchus cores, IRAS sources, and 
the Sco-Cen OB association. They find that the 
region contains some 23 13 CO clouds containing 
51 13 CO cores, with a total mass of 4OOOM0 and 
typical core densities of Nn 2 ~ 3 x 10 3 cm~ 3 . 
The dense cores (Nn 2 ~ 10 4 cm~ 3 ) and velocity 
str ucture were subsequent ly followed up in C 18 Q 
bvlTachihara et al.l (|2000aHbl . l2002l) . INozawa et ail 
(|199ll ) identify only thirteen YSOs associated with 
the cores, pointing to a low star formation effi- 
ciency of 0.3%. 

Distance estimates for the Oph N molecular 
cloud complex come from its relationship with the 
stars in Upper Sco (US) as, from extinction, the 
molecular clouds lie in fro nt of and distribute d 
through the OB population ( de Geus et al. 1989). 
Hipparcos parallaxes place Upper Sco at a mean 
distance of 145 ±2 pc ( de Zeeuw et al.|[l999l) . with 
a line-of-sight extent of ±17 pc assuming the 14° 
spatial extent is reproduced in the third dimen- 
sion. This places an effective upper limit on the 



clouds' distances of 162 pc. Extinction-based dis- 
tance modulus estimates suggest the clouds are 
distributed between 80 p c (near side) and 17 pc 
(far side), centre 125 pc ( de Geus et al.lll989h or, 
slightly fu rther away, 12 pc (near side) to 200 pc 
(far side) (|Straizvslll984 ). These distances are con- 
sistent with t he Hipparcos d a ta. L ooking at in- 
dividual stars, Ide Geus et al. (Il990h suggest that 
the western clouds (our OphN 4,5,6 their Com- 
plex 2) lie in front of \ Oph (150 pc) and the 
northeastern (OphN 1/2, Complex 4) in front of 
£ Oph (130 pc). Recent estimates for the better- 
studied L1688 Ophiuchus cl oud range from 120- 
145 pc (jWilking et al.ll2008l) . There is certainly 
no reason to believe that the Ophiuchus clouds 
all lie at the same distance, but for convenience 
we assume a working distance of 130 pc, which is 
consistent with the above estimates. 

In this paper, we present mid-infrared, Spitzer 
Space Telescope observations of the high column 
density regions of Ophiuchus North. The observa- 
tions and data reduction are described in Sect. [2] 
Results, including source statistics, YSO candi- 
dates, extinction maps and large-scale emission, 
are given in Sect. [3] with comments on individual 
regions in Appendix [A] The results are discussed 
in Sect. 2] and summarised in Sect. [SJ 

2. Description of Observations 

We observed Ophiuchus North in the mid- 
infrared as part of the Spitzer legacy program 
"Gould's Belt: star formation in the solar neigh- 
bourhood" (SGBS). The clouds were mapped un- 
der the working name 'Scorpius' (Sco) and ap- 
pear under this name in the SGBS catalogues. 
We present them here as the Ophiuchus North 
(OphN) clouds i n line with previous nomenclatu re 



( Nozawa et"ai1 Il99lt iTachihara et al.l l2000bllal . 



2002) which reflects their location predominantly 
within the constellation of Ophiuchus. Only our 
regions OphN 5,6 and LDN 43 lie beyond the 
constellation boundary in Scorpius. 

The SGBS program aimed to complete the 
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mapping of local star formation started by the 
Spitzer "From Molecula r Cores to P l anet-f o rmin ; 
Disks" (c2d) project (|Evans et al.l I2003L Ex} 



by targetting the regions IC5146, CrA, Scor- 
pius (renamed Ophiuchus North), Lupus II/V/VI, 
Auriga, Cepheus Flare, Aquila (including Ser- 
pens South), Musca, and Chameleon to the 
same sen sitivity and using the s a me reduction 
pipeline (iGutermuth et all l2008t lHarvev et al 



2008; lKirk et al.ll2009l:|Peterson et al.l201ll:ISpezzi et al 



20111 ). Images were made at 3.6/4. 5/ 5.8/8. 0/zm 



with the Infrared Array Camera (IRAC: lFazio et al 



20041 ) and 24,70 and 160/zm with the Multi- 
band Imaging P hotometer for Spitzer (MIPS; 
Rieke et al.ll2004l) . With an 85cm mirror, IRAC 
observes with an angular resolution of 2" whereas 
MIPS is diffraction limited with 6", 18" and 40" 
resolution at 24, 70 and 160/zm respectively. For 
our observations, we targetted small regions en- 
compassing the Ay > 3 contours from the op- 
tical extinction map of iDobashi et al.l (|2005l) . as 
shown in Fig. [TJ The area in Oph N which lies 
above Ay > 3 is fragmentary and scattered over 
an area of 20 sq. deg. Some of these regions 
(L158,L204,L146/CB68,L234E,L260,L43) had al- 
ready been mappe d as part of the S pitzer "Cores 
to Disks" project (jEvans et al. I l2003l) . These were 
avoided by the Spitzer Gould Belt team to avoid 
unnecessary duplication of observations, as the 
two projects work to the same target sensitivi- 
ties. Most of these c2d data are incorporated in 
this study of Oph N. The ex ception is L4 3 whic h 
is presented separately by IChen et al.1 (l2009h . 
Ultimately, seven new areas wer e mapped by 
SGBS with IRAC and MIPS. The IDobashi et all 
(|2005l ) and iRowles fc Froebrichl ([2009J) extinction 
maps, which are not limited to the IRAC ob- 
servations but extend across the entire area cov- 
ered by MIP S, confirm that all A y > 3 (mea- 
suring from Dobash i et al.l 120051 ) or Ay > 4 
( Rowles fc Froebrichl 120091 ) regions in these fila- 
ments were observed by IRAC with the exception 
of two small clouds 0.3° to the south of OphN 6 
(included in the MIPS map) and 0.4° to the north 
of CB68 (|Tachihara et al.ll2000bl core q2). De- 



tails of the datasets included for each region are 
listed in Table [TJ including the observation dates, 
Astronomical Observation Request (AOR) iden- 
tification numbers, program identification (30574 
for Spitzer Gould Belt, 139 for "Cores to Disks"), 



and duration. Table |2) gives the associated Lynds 
Dark Clouds JwdJn962) and molecular cores 
(C 18 0, iTachihara et al.ll2000al ) for each region. 

The overlap area covered in all four IRAC 
bands is slightly smaller than the area covered 
in any single band because the array design leads 
to an offset between the 3.6/5.8/im and 4.5/8.0^m 
maps. The final area covered by MIPS is much 
larger than with IRAC because of the long MIPS 
scans, with almost all the IRAC areas covered by 
MIPS as shown in Fig. Q] In total, 2.1 square de- 
grees were covered by IRAC and 6.5 with MIPS. 
This is roughly half the area of 14.4 sq.deg. cov- 
ered by MIPS in the main Ophiuchus clouds 
([Padgett et al.ll2008l) . 

The observations of each area were split be- 
tween two epochs to allow removal of transient ob- 
jects, with the second epoch maps offset from the 
first epoch maps so that bad pixels do not always 
lie in the same sky positions. The IRAC observa- 
tions were taken with a total integration time of 48 
seconds per point, split equally between the two 
epochs, and an offset of 10" in array coordinates 
between epochs. Short integrations in High Dy- 
namic Range (HDR) mode were also taken for all 
regions in which bright YSOs were expected (the 
exception was L234E). The SGBS MIPS observa- 
tions were taken in fast-scanning mode, stepping 
by 240" cross-scan to fill gaps in the coverage, and 
(125", 80") between the two epochs in order to pro- 
vide full 70/im coverage with only half the array 
working. MIPS total integration times were 32.4s 
at 24 and 70/im and 6.2 seconds at 160/zm. The 
"Cores to Disks" observations of the small cores 
L234E and CB68 were taken in MIPS photome- 
try mode. Each core was observed in 2 epochs. 
At 24/im 1 cycle of 3 seconds was taken at each 
epoch for a total integration time of 84s. At 70/im 
3 cycles of 3 seconds were taken at each epoch for 
a total integration time of ~ 94s. 

2.1. Data reduction 

Data reduction was carried out using the c2d 
pipel ine as described i n the d elivery documenta- 
tion dEvans et al.l 2007; see also lHarvev et al. 2006 
and iRebull et ail 120071 for details of IRAC and 
MIPS processing). The Basic Calibrated Data 
(BCD) from the Spitzer Science Center were pro- 
cessed by the SGBS team to remove artifacts ( 
eg. bad pixels, jailbar effects, muxbleed column 
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pulldown due to bright sources) and apply the 
location-dependent photometric corrections. 

2.1.1. Mosaics 

Mosaicking was carried out using the SSC 
"Mopex" code. For IRAC mosaics, the high dy- 
namic range (HDR) data were included in the 
final map, which improves the dynamic range and 
allows for the inclusion of otherwise saturated 
sources at the expense of slightly increased noise 
levels. For both IRAC and MIPS, both epochs 
were combined. 

The SGBS Oph N data were reduced to form 
six separate maps, numbered OphN 1-OphN 6 
in order of decreasing Galactic longitude. These 
were originally mapped as Sco 1-6 but relabelled 
'OphN' in line with the overall cloud nomencla- 
ture change. OphN 1,4,5 and 6 each incorporated 
a single SGBS AOR area. For OphN 2, the c2d 
cloud L204C was mosaicked with the SGBS data. 
For OphN 3, the MIPS observations for two of the 
SGBS regions and the overlapping c2d region L158 
were merged into a single mosaic. The IRAC ob- 
servations for OphN 3 remain split into OphN 3a 
(including L158) and OphN 3b. The two separate 
c2d cores L234E and CB68 were mosaicked indi- 
vidually. These pairings are indicated in Table [TJ 
A finding chart showing the six SGBS regions and 
the two compl e menta ry c2d areas overlaid on the 
iDobashi et al. I (|2005h extinction map is shown in 
Fig.Hl 

2.1.2. Catalogues 

Source extraction was carried out in each of 
the 4 IRAC bands on the combined mosaics (2 
epoch plus HDR), plus MIPS 24/im on the com- 
bined epoch mosaics, using c 2dphot, a photome- 



try tool based on Do PHOT (jrlarvev et all 120061: 



Schechter et al.lll993l ). Sources believed to be real 



at one or more wavelengths were band-filled at the 
missing wavelengths, first in IRAC bands 1-4 then 
MIPS 24/im then MIPS 70/mi by fitting a point 
spread function (PSF) to images at the missing 
wavelengths, and these fluxes are included in the 
catalogues with bandfilling noted in the data qual- 
ity flags. Fluxes were extracted from the 2MASS 



wavelength catalogue with an 8" position match- 
ing criterion. With a larger pixel size of 10" 
at 70/im (compared to 1.2" for IRAC and 2.55" 
at 24/itm) sources often matched more than one 
shorter-wavelength source. The 70/mi fluxes were 
assigned to the best candidate by hand matching 
the spectral energy distribution (SED), with the 
remaining possible 70/im emitters also noted in the 
catalogue. Source matching was not attempted 
at 160/im as the spatial resolution is low at this 
wavelength (~ 40") and most bright regions are 
saturated. For further details on source extraction 



see t he c2d delivery documentation (jEvans et al 
20071 ). 



point source catalogue (jSkrutskie et all 120061 ) us- 
ing a 2" position matching criterion. The 70/im 
source list was separately matched to the shorter 



2.1.3. 160 micron data 

The 160/im data obtained in fast scan mode 
mapping does not have enough redundancy to fill 
in completely all the gaps due to a dead readout 

and the i ntermediate gap between the array detec- 

tor ro ws ( MIPS Instrument and MIPS Instrument Support Teams! 
20071 ). Furthermore with only 2 epochs, effectively 
2 pointings per pixel, the effects of hard radiation 
hits and saturation translates into small regions 
without data. To deal with these gaps and to pre- 
serve as much as possible of the diffuse emission, 
the images are first resampled from 15" to 8" pix- 
els and then a 5 pixel by 5 pixel median filter is 
applied to the image. The net effect is a slight 
redistribution of surface brightness (~ 15%) and 
smearing of the original beam from 40" to about 
1' in size. 160/im data are not available for the 
two c2d cores CB68 and L234E. 

Fluxes at 160/im are determined using aper- 
ture photometry with a 32" radius aperture, an 
annulus from 64-128" and an aperture correc- 
tion of 1.97. Flux density uncertainties are about 
20% below 5 Jy and 3 0% for higher flux densi- 
ties (|Rebull et al.ll2010h . Visual inspection of each 
160/im source in the original unsmoothed image 
is used to determine whether the object is cleanly 
detected or contaminated by data dropouts. 

3. Results 

Red-green-blue (RGB) images for the eight 
Ophiuchus North regions mapped (OphN 1-6, 
CB68 and L234E) are shown in Figs. [land [3] The 
first figure shows a combination of IRAC 3.6/im 
(blue), 4.5/im (green) and 8/iin (red), and the 
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second combines IRAC 4.5/zm (blue), IRAC 8/im 
(green) and MIPS 24/zm (red). The RGB images 
each cover the overlap region where data is avail- 
able from all three chosen Spitzer instruments, 
corresponding to the IRAC areas (white boxes) in 
Fig. [T] The regions are irregular with sizes start- 
ing at 0.3 pc for CB68 and L234E, with OphN 3a 
the largest region mapped at a couple of parsecs 
in length; scale bars of 0.2 pc are shown on Figs. [2] 
and [3] 

Unextincted foreground stars appear blue- 
white in these images. Embedded protostars 
have redder colours, as do background galaxies, 
background giants, and planetary nebulae, ex- 
tincted ma in-sequence stars and foreground late- 
type stars ( Robitaille et al.ll2008l) . The red wisps 
of 8fim emission are PAH emission from heated 
clouds. Shocked H2 emission from molecular out- 
flows can show up at 4.5/xm in green on the IRAC 
3-colour plot, but there is no evidence for this 
here. 

Figs. [2] and [3] shows that there are no dense 
protostellar clusters in Ophuchus North, such as 
Spitz e r found in Serpens and Aur iga (jGutermuth et al 
20081: iMatthews fc Tothill|[2012h . The bright re- 
gion in the southwest of OphN 6 is the o ptical re- 
flection nebula IC4601 (|Magakianll2003l ). excited 
by a few late-type B and A stars (see Sect. IA.7I) . 
There is no dense young reddened cluster as- 
sociated with these intermediate-mass stars, as 
Fig. [3] shows. The young stellar objects associated 
with OphN 6 and other regions are discussed in 
Sect. [331 



3.1. Source statistics 

The Ophiuchus North IRAC and 2MASS source 
detection statistics are given in Table [3J In the 
area covered by all four IRAC bands, just over 
100,000 sources were detected in at least one IRAC 
band with a signal-to-noise of at least 3, with more 
than 5000 detected in all four IRAC bands. Nearly 
all of the 2MASS sources in the field were also 
detected by IRAC, with half of these detected in 
all 4 IRAC bands. The number of sources detected 
as a function of the required signal-to-noise in the 
individual wavelength maps is given in Table HI 
These counts include the entire region covered in 
each waveband and not solely the overlap area, 
which is why the 3.6/zm detection count for S/N=3 
is larger than in Table |3J 



3.2. Background sources 

Most of the sources detected by Spitzer are ex- 
pected to be either main-sequence stars or back- 
ground galaxies. In order to estimate the fraction 
of young sources associated locally with the Oph N 
star-forming region, a comparison with the Wain- 
scoat model for the Galactic stell ar background at 



this latitude is made in Fig. 2] (jWainscoat et al 



1992). The figure shows differential source counts 
for all sources in each of the 4 IRAC bands and 
MIPS1 from the Oph N catalogues (grey curves for 
all sources; black curves for sources classified as 
stars with 3cr detections) compared to the Wain- 
scoat predictions (dashed curve). For Oph N, most 
of the source detections are main-sequence stars, 
as shown by the lack of excess counts over the 
Wainscoat model in the IRAC bands. The ex- 
cess in the MIPS band over the Wainscoat model 
is mainly due to background galaxies. The com- 
pleteness limit at approximately 15 magnitudes in 
the IRAC bands can be seen in the turnover in the 
star counts (black line). 

3.3. YSO candidates 

The main aim of the SGBS programme is to 
locate the young stellar objects associated with 
star formation. In the Oph N molecular clouds, 
eleven candidate Young Stellar Objects (YSOc) 
were identified on the basis of their position in 
IRAC colour-magnitude diagrams, as shown for 
each of the IRAC 4-band detections in Fig. [5] 
The selection proc ess is described in detail in 
Harvev et al. (2007) but in summary, each colour- 
magnitude diagram yields a probability of a source 
being a YSOc with the region occupied by galax- 
ies, based on the SWIRE sample, giving a low 
probability. The individual probabilities are then 
combined with some additional criteria to separate 
the candidate YSOs (coloured red) from galaxies. 
The method aims to eliminate background sources 
at the expense of also eliminating some YSOcs. 
This process is not infallible and followup visual 
inspection of the mosaics identified one extended 
and elliptical YSO c andidate as a gala xy (6dFGS 
gJ164828.8 -141437: [Jones et al.ll2006h . 

The remaining ten YSOc were classified Class I, 
flat- spectrum, II, or III according to their spectral 
index fit between K-band (1.2/zm) and MIPS band 
1 (24/xm) as listed in Table [SJ Fluxes for these 
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sources are listed in Table [6] Spectral energy dis- 
tributions (SEDs) for these sources are shown in 
Figure [6l These include optical (B,V,Rc) mag- 
nitudes from the NO MAD catalogue from a 5" 
radius position search ( Zacharias et al.l|200 5f). For 
the more evolved sources, the SEDs are also shown 
dereddened to fit a stellar photosphere at the short 
wavelengths. The best wavelength/band for iden- 
tifying the continuum of stellar emission in these 
young objects is J-band, which was used to nor- 
malize the reddened photospheres to the SEDs. In 
all but two cases this was a K7 star, but for sources 
3 and 4 an AO star was needed in order to re- 
produce the high fluxes. These YSOc SEDs show 
disk excesses over a stellar photosphere at long 
wavelengths. The poor match between the SED 
and the optical data in some cases (YSOc 3,6,7) is 
most likely to be due to misassociations as several 
optical sources within the 5" radius contribute to 
the SED of a single source in the Spitzer bands. 
Some of these sources could also have ultraviolet 
excesses due to disk accretion that prevent good 
fits to main sequence stel lar photosphere models 
( Calvet fc Gullbrine 1998|) . It is also possible that 
the reddening law in these dense cores could differ 
from the Ry — 5.5 Weingartner b, Draine inter- 
stellar grain model (Weingartner fc Draind 12001 



Indebetouw et al ] l2005l) . but this correction would 
be likely to be similar for all sources. 

No photosphere is fitted to the flat-spectrum 
or Class I sources. These embedded YSOcs 
have emission peaking at longer wavelengths than 
main-sequence stars and can be identified by their 
redder colours in the RGB diagrams (Figs. [2] and 
[3D. The Class I sources in OphN 1 (YSOc 1) and 
CB68 (YSOc 10), and the flat-spectrum source in 
OphN 6 (YSOc 8), show up in this way. 

More than half of the YSOcs were detected 
by IRAS and all but three were previously iden- 
tified as YSOs (see Table p. The two deeply 
embedded Class 0/1 sources are the known pro- 



tostars L260 SMMld Andre fc Montmerld 



Visser et al.l|2002h and CB68 (|Carballo et al 



1994; 



1992; 



Huard et al.l ll999). Most of the Class II objects 



are T Tauri stars which have been the subject of 
multiple studies. The three new identifications 
are the Class III sources OphN YS06, 7 and 9, 
identified through their small 24^m excesses. Two 
of these (YSOc6 and 7) may be background AGB 
stars; we examine this possibility in Sect. 



OphN 7 and 9 both lie between the extinction 
peaks of OphN 6 (Fig. UTI) . Finding pre-main- 
sequence stars in OphN 6 is unsurprising as it lies 
on the edge of the Sco-Cen association. 

3.4. 24/im emission and MlPS-only YSO 
candidates 

A larger area of Oph N was covered by MIPS 
at 24, 70 and 160^m and by 2MASS in the near 
infrared, but not at intermediate wavelengths by 
IRAC. YSOs in that region are not found on the 
standard IRAC-based colour-magnitude diagrams 
as shown in Fig. [SJ A start on identifying YSO 
candidates with no IRAC data can be made by 
looking for sources with red colours in 2MASS 
K s - [24]. 

The Ks vs. Ks — [24] colour-m agnitude dia- 
gram is shown in F ig. [TJ following I Rebuff et al.l 
(|2007l) (Perseus) and lPadgett et all (|2008l ) fQphi- 
uchus). All Oph N sources detected at these wave- 
lengths are included (OphN 1-6, CB68 and L234E 
combined). The source symbols are based on the 
c2d pipeline classification. To be classified as a 
YSOc requires fluxes in all four IRAC bands and 
MIPS1; red sources which do not fulfil this crite- 
rion (eg. those not observed by IR AQ tend to 
be cl assified as 'gale' or 'star+dust' ([Evans et al 
20071 : see Sect. [2X2]). 



In Fig. [71 sources with the colours of stellar pho- 
tospheres appear to the left with Ks — [24] < 2, 
and most of the sources in this region are already 
classified as stars (grey circles) by SED fitting. A 
reddening vector is shown in the bottom left of the 
plot, assuming an extinct ion law of A2i/AK a = 
0.5 as recommended by IChapman et al. (2009); 
however, the 24/im relative extinction is quite un- 
certain. Reddened stars are likely to appear down 
and to the right of the main stellar population. 
The extinction in Oph N is low over most of the 
cloud, with only a small fraction above A v = 3 
(Figs. Q] and Sect. 13.71) . so few background stars 
should be reddened significantly, and indeed there 
are relatively few sources in this region. The lower 
left region of the diagram is not populated because 
stellar photospheres with Ks > 10 do not lie above 
the MIPS 24/im sensitivity limit unless they are 
reddened. 

Background galaxies ('Gale', green x) typi- 
cally have low Ks fluxes and lie towards the bot- 
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torn of the plot, as shown by the colours of ob- 
jects in the SWIRE catalogue ([Rebull et al.ll2007 : 
Padgett e t al. 2008). In this region are found most 
of the c2d pipeline-classified 'star+dust' sources 
(red +), which have the colours of a stellar pho- 
tosphere but excess due to dust at the long wave- 
lengths. 

Candidate YSOs with Kg — [24] excesses appear 
in the centre of the plot, to the right of the stars, 
depending on the contribution of dust in the disk 
and/or envelope. This is where the ten 'YSOc' 
already identified by the c2d pipeline and IRAC 
colour-magnitude diagrams appear (red filled cir- 
cles), with the exception of CB68, which was not 
detected by 2MASS (or IRAC1) and does not ap- 
pear on this plot. In this region there are several 
additional YSO candidates with no IRAC detec- 
tion, currently classified 'star+dust' (red crosses) 
due to their long-wavelength excess over a stellar 
photosphere. 

At this point we make a selection for a clean 
sample of MlPS-only YSO candidates, reject- 
ing the regions of colour-magnitude space with 
Ks < 13, which excludes the majority of back- 
ground galaxies, Ks — [24] > 2, which excludes 
the background stars, and [24] > 9 which re- 
moves reddened stars. This is a conservative 
selection which probably rejects some genuine 
YSOs. We make the gala xy cut at Ks = 
13 rather than K s = 14 dRebull et al.l 12007; 



Padgett et al.1 I2008D because the two brightest 
sources below Ks — 13 in the flat-spectrum 
region of the plot both have galaxy identifica- 
tions. The misclassified 'YSOc', which has an 
IRAC detection and appears as a red filled circle, 
is 6dFGS gJ164828.8-141437 and the brighter 
K s ~ 13 source is 6dFGS gJ163139.8-202044. 

The five MlPS-only YSO candidates in Oph N 
selected by these criteria are listed in Table . 
We can classify them further based on the stan- 
dard spectral index a calculated from Ks — [24] 
into Clas s I, flat-spectrum, Class II and Class II I 
sources (|Rebull et al.l 120071: iGreene et ail Il994[) ; 
the corresponding regions of colour space are iden- 
tified on Figure [7] The three Class II candidates 
M5, M2 and M4, in OphN 3,4 and 6 respectively, 
are close in colour to the existing YSOcs and have 
strong Ks and 70/zm detections. The fainter Ks 
source M3 is also a Class II candidate in OphN 6. 
In OphN 3, the Class III source Ml which lies red- 



wards of the main stellar population is identified 
as IRAS 16484-1557. 

To fill in the missing mid-IR portion of the spec- 
trum for regions not covered by IRAC, additional 
mid-infrared fluxes were taken from the WISE 



Preliminary Release Source Catalog ([Wright et al 
2010l fl. WISE provides fluxes at 3.4, 4.6, 12 



and 22/Ltm with angular resolution of 6-12". The 
resulting spectral energy distributions, shown in 
Fig. [51 confirm the sources to be similar to the 
Class 11 /III YSOcs in Fig. H The classification 
of the source M5 (V1121 Oph), which shows a 
strong disk excess, deserves some further com- 
ment. Aperture photometry with a 100" aper- 
ture gives a 24/mi flux an order of magnitude 
higher than the c2d pipeline (~ 4000 rather ~ 
400 mjy) and the higher flux is confirmed by 
WISE at 22/im (5915 ± 49 mjy); see Table. [8]). 
The higher 22/24/zm fluxes move this source to 
the Class II/flat-spectrum boundary (Fig. [7]). 

Ks — [24] alone is an imprecise selection tool, 
and there may be additional YSOs in the ex- 
cluded regions, particularly among the faint 24/im 
sources. Colour-magnitude searches based on 
WISE data may be a way to identify these, but 
are beyond the scope of this paper. 

3.5. AGB star contamination 

There remains the possibility that some of 
the YSOcs may be background asymptotic gi- 
ant branch (AGB) stars. Like YSOs, AGB stars 
have red colours due to the dust shells which sur- 
round them. They occupy a similar region in 
colour-magnitude space to YSOcs, and cannot 
easily be separated by the criteria shown in Fig. [5] 



(Robitaill e et al.l 12008). As befits an old stellar 



population, AGB stars are distributed through- 
out the Galaxy with increased counts in the thick 
disk and bulge. Oph N lies 17-24 deg from the 
plan e, so using the Galac tic distribution model 
from llshihara et all ( 2011 ) an average AGB con- 
tamination of less than 0.5 deg -2 would be ex- 
pected, or fewer than 1 AGB star in the 2.1 deg 2 
area covered by IRAC and fewer than 3 AGB stars 
in the 6.5 deg 2 regions covered by MIPS. 

AGB stars can be identified by their low proper 
motions, reflecting their Galactic distances and in- 
consistent with membership of nearby Upper Sco. 

1 http: / / wise2.ipac.caltech.edu / docs / release/prelim/expsup / 



8 



From Hipparcos, Upper Sco has average proper 
motions o f (ua,us) = f~ll'~2 4 ') in Galactic co - 
ordinates ( de Zeeuw et al. 1999t Mamaiekl [2008) . 
Two of our IRAC+MIPS1 selected YSOc have 
inconsistent and low proper motions: YSOc 6 
(IRAS 16191-1936 in OphN 3b) ha s proper mo 



tions of 4.7 ± 7.2, -1.2 ± 7.2 mas/yr (|Roeser et al 



2010h and YSOc 7 (OphN 6) has proper motions 
of —2.6 ± 5.5, 2.3 ± 5.5 mas/yr. Neither of these 
are consistent with a nearby location in Oph N, 
whereas the other YSOcs all show proper motions 
consistent with Upper Sco. 

Additionally, AGB show a steep spectral in- 
dex at long wavelengths: a criterion of [8.0] — 
[24] < 2.2 can be used to select them, though 
this criterio n alone will not distinguish them 
from YSOcs (|Whitnev et al.ll2008l : lRobitaille et al 
2008). Fig. [9] shows several combinations of 



colours and magnitudes of the YSO and AGB can- 
didates compared to those of the known AGB sam- 
ple in Serpens (iden tified by infrared spectroscopy; 
Oliveira et aj]|2009h . It can be seen that the two 
IRAC-observed AGB candidates in Oph N lie in 
the same region in colour-magnitude space as the 
confirmed AGB stars. YSOc 6 (IRAS 16191-1936 
in OphN 3b) has [8.0] - [24] = 1.2 and YSOc 7 
(OphN 6) has [8.0] - [24] = 1.0. These two sources 
are among the most luminous in our sample, as 
shown in Fig. El In Fig. [9j YSOc 9 also has 
the colours of an AGB candidate, but its proper 
motions are consistent with membership of Up- 
per Sco din-, Us = -3.9 ± 4.9, -29.9 ± 4.9 mas/yr; 
Roeser et alll2010h . 

In addition, two of the MlPS-selected YSO can- 
didates are likely to be background AGB stars. 
Sources Ml and M2 both have low proper mo- 
tions inconsistent with Upper Sco: Ml (fi a ,Us = 
-4.7 ± 4.9, -1.0 ± 4.9 mas/yr) and M2 (-4.7 ± 
4.9. -1.0 ±4.9 mas/vn lRoeser et alJ l2oToh. IRAC 
fluxes at 8/im are not available to apply the Whit- 
ney et al. [8.0] — [24] < 2.2 colour criterion exactly, 
but applying the same principle to WISE fluxes 
these two sources have the steepest spectra of the 
sample between 12 and 22/im ([12] - [22] < 1.0), 
supporting their identification as AGB stars. In 
the ASAS variable stars catalogue, Ml is identi- 
fied as a Mira variable with a magnitude varia- 
tion of 1.9 mag and period 274.050380 days (ASAS 
165122-1602.9. IPoimanski et al.l[2005h . 

The four AGB candidates among the YSOcs 



are identified by footnotes in Tables [5] and [7] and 
in Figs [T0HT3l They are not counted in our final 
YSOc tally (Sect.g). 

3.6. 70 micron emission and SCUBA maps 

Nine out of the ten YSOcs listed in Tableland 
two out of four of the MlPS-only YSOcs in Table[7] 
were observed and detected at 70//m. Detections 
at 70/im are either a sign of disk excess or a dusty 
envelope, as is clear from the SEDs (Fig. [6]). Often, 
these sources can also be detected in the millime- 
tre/submillimetre, on the long- wavelength side of 
the peak of the SED. In Fig.^Ul we show the MIPS 
70/im maps overlaid with contours of submillimc- 
tre emission where available. The 70/im maps 
show stri ping artifacts rema ining from the recon- 
struction ([Evans et al. 2007 ). We searched the re- 
reduced SCUBA archive for 8 50£im maps associ- 



ated with the Oph N regions (|Di Francesco et al 



2008) and found existing small SCUBA maps 
associated with parts of OphN 1 and OphN 3 



originally publis hed in lVisser et al.ll2002D. L234E 



originally published m lvisser et 
..Kirk et al.ll2005l) and CB68 jHuard et al. 
Vallee et al.l 120031: I Young et al.l l2006h . The two 



19991: 



YSOcs in the south of OphN 3 and in OphN 6 
do not appear to have been mapped by SCUBA, 
and we are not aware of any other mm/submm 
continuum observations of these sources. 

Of the two Class I YSOcs, CB68 has a 70/im 
detection: the OphN 1 Class I lies just off the edge 
of the 70/im map. Both have been mapped and 
detected at SbOum by SCUBA (jVisser et al.ll2002 : 
iHuard et al. 1999h . The flat-spectrum source 
(YSOc 8 in OphN 6) is also detected at 70/im. 
All but one of the Class lis are detected: the ex- 
ception is the MlPS-only Class II YSOc M4 in 
OphN 6. Neither the Class III source YSOc9 nor 
the AGB candidates YSOc7 and Ml are detected 
at 70/im, as expected from the mid-infrared falloff 
of the SEDs (Fig. ©and©. 

3.7. Extinction maps 

Extinction maps provide an alternative mea- 
sure of the cloud structure in Oph N. The Spitzer 
Gould Belt catalogue includes a measure of the 
visual extinction A v towards every source with a 
SED fitted by a reddened stellar photosphere (i.e. 
classified as a star) . By convolving the extinctions 
at each position with a Gaussian beam, these val- 
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ues have been turned into extinction maps for all 
the area in Ophiuchus North covered by IRAC. 
The detailed proc edure is explained in the c2d de- 
livery document ( Evans et al.l 120071 ) . but in sum- 
mary, extinction maps are provided at a range of 
spatial resolutions, with the maximum resolution 
limited by the number of stars which contribute 
extinction measurements within each beam. As 
this number decreases, the uncertainties in the 
extinction increase. For the Gould Belt extinc- 
tion maps to be accepted, a maximum of 10% 
of beams at any given A v level were allowed to 
be undersampled, defined as containing no stars 
within the beam FWHM. Thus the lowest resolu- 
tion Oph N maps were sampled with a 300" beam 
but the highest resolution A v maps vary from 90" 
(OphN 4 and 5) through 120" (OphN 1,3 and 6) 
to 150" in the high-extinction OphN 2. The ex- 
tinction maps are made primarily using stars with 
IRAC detections and at least one good 2MASS 
detection. In addition, a limited number of stars 
with no 2MASS detection but good IRAC detec- 
tions are used to supplement the measurements in 
the higher extinction regions (Ay > 8.5), where 
the number of near-infrared detections is reduced. 

Extinction maps at 150" resolution (the high- 
est resolution available for Oph N regions 1-6) 
and 120" resolution (L234E, CB68) are plotted 
as black contours over the MIPS 24/im mosaics 
in Fig. [TT] The advantage of the Spitzer-derived 
extin ction maps over the optical (iDobashi et al 



20051 ) or near-infrared (|Rowles fc Froebrichl 120091) 
is that they can probe high extinctions at rel- 
atively high resolution and so can detect com- 
pact, high column density cores. This is particu- 
larly useful in Oph N where extensive mm/submm 
maps are not yet available. 

The mid-infrared-derived extinctions shown 
in Fig. [TTJ are higher than the extinctions de- 
rived from 2MASS colours ( Rowles fe Froebrichl 
20091 : iFroebrich et al.l 120051) by typically 2 mag- 
nit udes and from optic al star-count extinctions 
of Dobashi et al. (|2005l ) by typically 3 magni- 
tudes. This can be seen from the contours on 
the extinction maps in Fig. [TT1 th e Ay > 3 re- 



gions from the lDobashi et al" 



20051 ) star-count ex- 



tinction maps (blue dashed c ontours) correspond 
rough ly to Ay > 4 on the iRowles fe Froebrich 
(2009) near-IR extinction maps (shown as red 
dashed contours) or Ay > 6 on the Spitzer ex- 



tinction map (base level of black contours). The 
IRAC selectio n (bla ck boxes) was based on the 
Dobashi et al.l (|2005l ) Ay = 3 contour. The rea- 



son for this difference in extinction level is not 
simple and requires further explanation. A sig- 
nificant difference between the extinction maps 
produced for the Spitzer Gould Belt survey com- 
pared to those for c2d is that no extinction offset 
is subtracted from the Spitzer Gould Belt regions. 
For the c2d maps, an extinction offset was calcu- 
lated using off-cloud fields believed to be free of 
extinction from the molecular cloud. These off- 
sets, which for the Ophiuchus cloud lying nearest 
to Oph N is Ay — 1.96, were subtracted from 
the on-cloud extinctions. It was believed that 
these corrections were for foreground, line-of sight 
extinction from the diffuse interstellar medium. 
From a detailed comparison of t he mid-infrared 



extin ction law with the models (|Chapman et al 



2009) it is now thought instead that these ex- 



tinction offsets are due to deviations from the 

5.5 



assumed IWeingartner fc Draind (|200ll ) Ry 
extinction law in the mid-infrared region of the 
spectrum when Ay is low. Ay was calculated as- 
suming the extinction law can be parameterised as 
A\(Ry,X) with Ry = 5.5 (IWeingartner fc Draine 



2001; Cardelli et al 



19891 ). The conversion of 
mid-IR extinction to Ay is nonlinear with the 
Ry = 5.5 model producing reasonably reliable es- 
timates of Ay at high extinctions (Ay > 10) but 
overestimates at low extinctions (Ay < 10) by a 
factor which varies with Ay and cloud. This is the 
cause of the apparent extinction in the off-cloud 
fields, and the difference of Ay = 2-3 from the 
extinction maps derived at shorter wavelengths. 

The high extinction regions are small and frag- 
mented, typically less than a parsec in size. The 
masses in each region above extinction thresholds 
of Ay = 6, 10 and 20, derived from these maps, 
are given in Table [9] Approximately 100 solar 
masses of gas lies at Ay > 10 in Oph N, split be- 
tween 12 clumps (Fig. [TT]). This is about 1/30 of 
the ma ss in Ophiuchus though similar to Lupus I 
or IV ( Heiderman et al. 2010l ). Each individual 
clump only has a mass reservoir of order 10M Q 
sufficient to form a few low-mass stars,. Several 
Oph N regions show cores of Ay > 20 in the 150"- 
resolution Spitzer-derived extinction maps: there 
is one in OphN 1, two in OphN 2, and three in 
OphN 3 (Fig. [TT]). The OphN 1 and CB68 cores 
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contain the two Class I YSOs OphN YSOl and 
10, respectively, but the other extinction peaks 
are currently starless. Class II and III YSOs lie 
on the edges of the extinction regions in OphN 3a, 
3b and OphN 6. 

3.8. 160/im and IRAS maps 

Emission at wavelengths of 100/^m and longer is 
dominated by cold dust and extended cloud struc- 
ture. In addition to the MIPS 160/mi data, we also 
obtained IRAS 100 jam images that were resolu- 
tion enhanced using the HIRES algorithm^. Im- 
age coordinates were additionally precessed from 
B1950 to J2000 coordinates using the Goddard 
IDL astron library. The algorithm improves the 
spatial resolution from the n ative IRAS (3.8' x5.4') 
to 2' at 100 fxm wavelength (|Cao et al.lll997l) . The 
MIPS 160/im data have a factor 5 better spatial 
resolution of 40". The HIRES processing pro- 
vides sharper and more clearly defined filaments 
and other extended structure. Point sources are 
however surrounded by a ringing artifact. Since 
our focus is on the extended structure the point 
sources are removed from the 100 /im image using 
median pixel replacement. 

The Spitzer 160/xm and IRAS 100/im maps are 
shown side-by-side in Fig. HU The IRAS 100/im 
maps have broader spatial coverage and show 
the continuous filamentary structure linking the 
OphN 1 to 3 clouds. The OphN 4-5 clouds are also 
linked by 100/im emission which takes the form of 
several cavities or bubbles. Although there are 
ionising stars in the region, these lie to the north 
(the B2V star X Oph at I, b = 357.93, +20.68 and, 
in t he past Myr, ( Oph ; see Fig. Q] of this paper 
and lNozawa et ah 1991 Fig. 5) whereas the open- 
ings of the bubbles lie to the south. The bubbles 
themselves contain no ionising stars to produce 
HII regions. A possibility is that they could be 
fluid instabilities at a hot gas / cold cloud inter- 
face where the flow from the Upp er Sco HI bubb le 
has broken out (see discussion in lde Geuslll992l ). 

The extended emission exhibits structure that 
is very similar at both 100 and 160 /im, apart 
from resolution effects. This similarity indicates 
that the two wavelengths are mostly tracing the 
same dust component. To determine the dust tem- 
perature and derive physical parameters we follow 



2 http: //irsa. ipac . caltech. edu/IRASdocs/hires_over .html 



the technique presented in iTerebev et al.l ((2009) 
and briefly described here. For thermal dust emis- 
sion the intensity I v at frequency v is given by 
I„ = t v B v (T), appropriate for low optical depth 
and where B V (T) is the Planck function and t v the 
optical depth. The wavelength dependence of the 
optical depth is given by the usual t v ~ v$ where 
(3 = 1.5 — 2 at long wavelengths. If, for example, 
the dust temperature is constant across the image 
then the Planck term B U {T) is constant. The re- 
sult is that the 160 and 100 /im images will look 
identical, apart from resolution effects, while the 
structure in the image will linearly trace the opti- 
cal depth i.e. column density. Also, a plot of 1100 
versus 1160 intensity values for each image pixel 
will exhibit a linear trend whose slope is related 
to the dust temperature. In practice the assump- 
tion of constant dust temperature works well over 
scales of 1 to 2 degrees. 

Table [TU] shows the 1100 versus 1160 slope and 
corresponding dust temperature for each of the 
seven OphN regions separately. There is a small 
but real trend in dust temperature across the re- 
gion. The coldest dust temperatures are found in 
OphN 1 (15.6 K) and OphN 3S (15.9 K), located 
on the eastern side, while the warmest tempera- 
tures are in OphN 5 (17.0 K) and OphN 6 (16.8 
K) located on the western side, nearest the OB 
association as seen in projection. The higher dust 
temperatures on the western side support the idea 
that the OB stars enhance the local interstellar 
radiation field. In comparison, the Taurus star- 
forming region has colder 14.2 K dust (Terebey et 
al 2009, Flagey et al 2009), consistent with the 
lack of luminous heating sources in Taurus. 

All but one of the Oph N regions shows evidence 
for cold clumps in the excess map (Fig. [T3| . In 
OphN 1, these are associated with L260 SMM1 
and 2. The 160/xm excess peaks at the Class 1 
protostar OphN YSOcl which is indicative of a 
cold disk or envelope. By contrast, there is no 
enhancement at the positions of the Class II and 
III sources in OphN 3 although the clouds in the 
south of the region appear cold. The effects of 
UV heating, presumably from £ Oph, can be seen 
along the northern edges of these clouds as bright 
100/im features and negative values in the 160/im 
excess maps. OphN 4 exhibits one of the series 
of bright shells or loops best seen in the 100/im 
maps. The moderate extinction clump on the edge 
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of this cavity is warm. There is also little cold 
dust in OphN 5. By contrast, OphN 6 shows two 
cold clumps and cold dust associated with the flat- 
spectrum source OphN YS08 (L1719B). 

4. Discussion 

The Spitzer data advance our knowledge of 
Oph N with an improved list of candidate YSOs, 
higher-resolution extinction maps, and dust tem- 
perature estimates. In this discussion, we apply 
this information to revisit the YSO population, 
the current star formation rate / efficiency, and the 
evidence for triggering and other effects of the 
neighbouring Upper Sco OB association on the re- 
gion. 

4.1. Star formation count and efficiency 

How many stars are currently forming in the 
Oph N clouds? In the areas covered by both 
Spitzer IRAC+MIPS, ten YSO candidates are 
identified. In the areas with no or limited IRAC 
coverage, MIPS+2MASS colours add a further 
five. A further Class I so urce, RN091, is detected 
in L43 ( Chen et al.ll2009l) . bringing the total to 16. 
This count a lmost doubles the nin e sources known 
from IRAS (jCarballo et al.lll992l) . However, four 
of these candidates (including two of the IRAS de- 
tections) are probable AGB stars based on their 
colours and proper motions (Sect. I3.5[) . Exclud- 
ing the AGB stars, the total number of Spitzer- 
identified YSOcs in Oph N is 12. 

However, the area of Oph N molecular clouds 
surveyed by Spitzer is limited and it is possible 
that additional YSOs lie outside the mapped re- 
gion. Based on the low extinction these are likely 
to be Class II/III sources. We can estimate the 
number of YSOcs missed by Spitzer by considering 
the number of IRAS detections within and outside 
the Spitzer regions. IRAS covers the whole of the 
Oph N region but is less sensitive than Spitzer. of 
the 12 Spitzer YSOcs (including RN091 in L43, 
excluding probable AGN), only seven have IRAS 
identifications. 

A wid er area search f o r YSO s in Oph N was 
made by ICarballo etahl (jl992l) based on IRAS 
colours. Excluding the Ophiuchus and Lupus 
clouds, they identify six IRAS sources in the 
Oph N region as definite YSOs, five of which 
are cross-identified by Spitzer and one of which 



lies outside the Spitzer regions (T Tauri stars 
V1003 Oph or RNO90 near L43). A further 
Spitzer YSOc (Oph N YSOc3 or IRAS 16455-1405) 
is among a further ten IRAS sour ces ambiguously 
classified bv lCarballo et al. (1992) as either YSOs, 
galaxies or EGOs and one further Spitzer can- 
didate, YSOc 4, has an IRAS identification not 
list ed by Carballo et al. (1992) (though included 
bv iNozawa et all 1991 ). bringing the total num- 
ber of Spitzer sources detected by IRAS to seven. 
Hence, with the count of IRAS YSOs outside the 
Spitzer regions numbering one, our Spitzer survey 
encompasses 7/8 or 88% of the IRAS-identified 
YSOs in the region. 

In a ddition, there are fou r ambigous (as classi- 
fied bv lCarballo et al.lll992h IRAS detections out- 
side the Spitzer regions which have associated ex- 
tinction Ay > 1. If these also turn out to be YSOs 
(eg. IRAS 16534-1557 or IRAS 16439-0945; see 
sects. IA~T1 and |A~8"|) the completeness decreases fur- 
ther to 7/12 or 58%. None of these are associated 
with MlPS-only sources. Assuming these statis- 
tics hold for all 12 of the Spitzer identified YSOcs, 
and not only the seven detected by IRAS, the true 
number of YSOs in Oph N (north of the Ophiuchus 
L1688 and L1689 clouds) is likely to lie between 
14 and 21. 

The star formation surface density in Oph N 
has b een addressed recently by iHeiderman et al 



(|2010l ) (under the region name 'Scorpius'), who 
find a star formation rate per unit area Ssfr = 
0.343 ± 0.20 M Q yr^kpc" 2 based on the ten 
Spitzer IRAC+MIPS identifications reported here 
and comparing with an Ay > 2 area of 7.29pc 2 . 
This is among the lowest in Gould belt clouds 
(mean Esfr = 1.2 ± 0.2 M© yr _1 kpc~ 2 ), compa- 
rable to Lupus I, Auriga N and IC5146 E/NW 
but only 15% of that in Ophiuchus (2.45 ± 
1.6 M Q yr^kpc -2 ). The inclusion of RN091 and 
the three MIPS+WISE identifications, but the 
exclusion of two AGB stars, increases the rate to 
Esfr = °- 412 ± °- 25 M Q yr _1 kpc" 2 . 

An estimate of the instantaneous star forma- 
tion efficiency (SFE) Af»/(M* + M gas ) can be 
made comparing the total twelve YSOcs to the 
total mass in the Oph N clouds, w hich is 3000M W 



calculated for Ay > 1 in the Dobashi et al 



2005 ) map, consiste nt with 4400 M© from 13 CO 



( Nozawa et al.Ml991 



an average 



Assuming 12 YSOcs with 
stellar mass of 0.5 M© (following 
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Evans et alj|2009h gives an overall SFE of 0.20%. 
Taking into account the survey incompleteness 
gives an upper limit of 0.34%. These efficiencies 
are consistent with t he 0.3% upper limit given by 
Nozawa et al. (|l99lh (we assume half their aver- 
age stellar mass and a slightly different star count, 
12 rather than 13). 

Nozawa et al.l (1991) suggest two reasons for 



the low SFE: the cloud is clumpy and fragmented, 
with only a small fraction of the mass in high 
column density cores; and the UV radiation es- 
timated from the nearby OB stars is a factor of 
2 higher than the standard Habing field, lead- 
ing to high cloud ionisation and magnetic sup- 
port. These estimates of the instantaneous SFE 
do not take into account the earlier formation of 
the main-sequence stars in Upper Sco . How ever, 
we confirm the result of Nozawa et al 1 (|l99lh that 
the current star formation efficiency in the cloud 
complex is very low. 

4.2. Comparison with Ophiuchus 

This low star formation efficiency is com- 
pletely unlike the nearby Ophiuchus L1688 cloud, 
which is forming a dense embedded cluster (e.g. 
Padgett et al. 20081 ). Measured in low-excitation 
13 CO 1-0, the Oph N and Ophiuchus cloud 
masses are similar: 4000M Q compared to 30 50M^ 
(L16 8 8, L1689 and L1712 -1729 filament) (jLoren 
19891 : Nozawa et al. 1991 ). The main difference 
between the c loud complexes is the mass in dense 
gas. From the lRowles fc Froebrichl l|2009l ) 2MASS- 
based extinction maps, calculated identically for 
the two clouds, we calculate that Oph N has 
25% the mass of Oph above A v = 3 (680M© 
in Oph N), dropping to 4% of the mass at Ay > 5 
and les s than 0.5% for Ay=lQ. iTachihara et al.l 
(|2000bh found that traced in C ls O 1-0, the cores 
in L1688 are ten times denser than the others 
in the region, and 32% of the cloud mass is 
traced by C 18 compared to 8% in the Oph N 
clouds, i ndicating substantial op tical depths in 



the 13 CO(lTachihara et al.1 2000b). As studies ex- 



tend beyond the more obvious star-forming re- 
gions, many molecular clouds are seen to fit this 
pattern of extended low column density material 
with a few low-mass cores, but the Oph N clouds 
are particularly inefficient: the Lupus molecular 
clouds - also on the edge of Sco-Cen association 
- have a star formation efficiency of a few percent 



(jTachihara et alJll996t Ikerin et al.ll2008l ): Taurus 
is also less extreme than O ph N with a star for - 
mation efficiency of 1-2% ()Mizuno et al. 1995). 
Maddalena's cloud (G216-2.5) is a more extreme 
example with a mass of more than a million Mq 
only a few dozen YSOs (jMegeath et al.ll2009h . 



4.3. Prestellar cores 

There is no evidence that star formation is com- 
ing to an end in the Oph N region. The YSO 
candidates cover the full range of evolutionary 
classes (3 Class 0/1; 1 flat-spectrum; 7 Class II; 
1 Class III), and there are also several dense star- 
less cores in the region, including submm cores 
in L234E, L255 SMM2 in OphN 1, and L158 in 
OphN 3a. Star formation is clearly still ongoing 
in Oph N, though at a low rate. 

A statistical estimate of the starless core life- 
time can be made from the high extinction cores. 
Considering Ay > 10 cores, the ratio of star- 
less to protostellar cores is 12:3 (from FigQTJ in- 
cluding CB68, L234E and one starless and one 
protostellar core in L43). If the starless cores 
are prestellar cores, then this ratio suggests their 
lifetime is a factor of 4 long er than the proto- 



stellar lifetime of .54 Myr (|Evans et al.1 12009L 
Hatchell et~aT1l2007l ) or - 2 Myr. This is consis- 



tent with_exjjectelio^^ volume den- 

sity ((Ward-Thompson et al.ll2007i Fig.2) given a 



typic al core density of 7 x 10 3 cm 3 (jNozawa et al 



1991L from C ls O). This is further support for the 
idea that the restrictive step in Oph N is the for- 
mation of dense cores from low density material. 

4.4. Influence of Upper Sco OB associa- 
tion 

With a known age gradient through the OB 
associations from Upper Centaurus-Lupus to Up- 
per Scorpius to Ophiuchus, Ophiuchus North 
would seem an ideal region in which to look for the 
effect of triggering on star formation. But while 
there is ample evidence for the dynamical and ra- 
diative influence of Upper Sco on the region, the 
low star formation efficiency and normal fraction 
of starless cores do not suggest that star forma- 
tion has been promoted, even in the clouds most 
strongly affected by the nearby OB association. 

An expanding HI shell surrounds Upper Sco, 
created by the stellar winds and possibly the su- 
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pernova explosion of the binary partner of run- 
away £ Oph, impacting on the Ophiuchu s and 
Oph N clo uds from t he we st (|de Geuslll992h . Ac- 
cording to Ide Geusl (|1992l) , the expansion of the 
shell in this region is halted by the Ophiuchus 
and OphN 6 dense cores, with molecular gas swept 
from the main clump to produce OphN 4-5 and the 
elongated Ophiuchus filaments to the south. A rel- 
atively recent origin is consistent with the absence 
of star formation in OphN 4 and 5 but could be 
tested by chemical comparison with OphN 6. Both 
OphN 5 and OphN 6 have elevated temperatures 
(Sect. I3.8j) . The bubbles seen in the 100/mi maps 
of OphN 4 and 5 also suggest dynamic disruption 
of these clouds. A similar though less extreme lack 
of star formation is seen in Oph L1689 compared 
to L1688, wh ere the effects of t he O-star a Sco 
are strongest ((Nutter et al.ll2006t l. The Oph N 6 
and Ophiuchus clouds are at similar distances to 
Upper Sco, so the physical input from winds and 
radiation must have been similar in each region. 
The difference in the star formation points to dif- 
ferences in the initial conditions, with the gas in 
Oph L1688 sufficiently gravitationally bound to 
survive - and even benefit from - the onslaught 
from Upper Sco, whereas the gas in OphN 6 had 
too low a density to counter the feedback forces. 
Timescales are too short for the formation of the 
Class II and III sources in the L43, OphN 3 and 
6 clouds to have been initiated by this most re- 
cent shock, which has an upper limit on its expan- 
sion timescale of 2.5 Myr and may not even ye t 
have reached these clouds ( Tachihara et al. 1996). 
There is also no clear sequence in protostellar ages 
from the front to the back of the filament which 
would suggest that the star formation has been 
triggered by the interactions with Upper Sco OB 
association. 

The UV field from the runaway O-star ( Oph 
forming the S27 HI I region is possibly an ear- 
lier i nfluence on the norther n Oph N clouds and 
L43 (jTachihara et al] l2000al) . The dense gas in 



OphN 1-2 mainly lies along the west side of the 
filament and the clouds are also heated along this 
side (Sect. I3.8|) . The 100/im maps show bright 
rims to the west and north in the direction of 
£ Oph and sharp boundaries in the dust distri- 
bution. ( Oph (09.5) originates in Upper Sco 
and has since tracked from southwest to northeast 
across the northern Oph N region on timescales 



of a few x 10 5 years (see iTachihara et al.l 12000a , 
Fig. 5). The radiation has impacted OphN 1 and 
L43 and may have contributed to the compression 
of the cores hosting Class 0/1 sources in L43 and 
L260. The 2 Myr old T Tauri stars in OphN 3 and 
RNO90 near L43 predate the impact of the HII 
region. A study of the energetics suggests that 
UV photodissociation from ( Oph is currently ac- 
celerating gas fr om the remaining dense cores to- 



wards the east ( Tachihara et al. 2000ah . This is 



also suggested by the bright rims and cometary 
nature of OphN 1-3 as seen in the 100/zm IRAS 
map (Fig. [T2|). It seems likely that the small star- 
forming molecular cores in Oph N are the surviv- 
ing remnants of the giant molecular cloud which 
produced Upper Sco and are now slowly being ab- 
lated by the massive stars. 

5. Conclusions 

We have surveyed the high extinction regions 
(Ay > 3) of the Oph N molecular cloud complex 
in the mid-infrared with Spitzer IRAC and MIPS. 

There is little active star formation in Oph N. 
Twelve YSOcs are identified in total. Eight 
of these are identified from their IRAC-MIPS 
colours. An additional Class I protostar exists 
in the L43 cloud in the centre of this region 



((Chen et al.l 120091 ). Three further YSO candi- 
dates are found based on red [K — 24] colours. 
Four sources initially selected as YSO candidates 
were rejected as AGB stars and one as a galaxy. 
Of the twelve remaining YSOcs, three are Class I, 
one flat-spectrum, seven Class II and one Class III. 
All except one Class III and two Class II sources 
were previously known. The region as a whole has 
a low star formation efficiency of < 0.34%. 

The high extinction regions in Oph N are frag- 
mented into twelve small (~ 0.2 pc), scattered 
and low-mass (~ 10M Q or less) cores, of the 
kind which might form one or two low-mass stars. 
Three of these (OphN 1, CB68 and L43) currently 
contain Class I sources. The remainder are star- 
less. 

The interstellar medium in this region is influ- 
enced by the Upper Sco subgroup of the Sco-Cen 
OB association. There is evidence for dynamical 
interaction in the OphN 4 and OphN 5 bubbles, 
elevated temperatures in OphN 5 and 6, and irra- 
diated cloud edges throughout the region. The 
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bulk of the gas mass is at low column density 
with only a few low-mass cores surviving to form 
a few YSOs. This is very different from the situ- 
ation in nearby Op hiuchus L1688, whic h contains 
hundreds of YSOs (IPadgett et alj|2008ft. but sim- 
ilar to the Lupus clouds (jTachihara et al.l 1 19961: 



Merin et al.l |2008) . As the impact of radiation 



and winds from Upper Sco must have been sim- 
ilar in both regions, the initial conditions must 
have been very different. Whereas the main Ophi- 
uchus complex contained hundreds of solar masses 
of dense gas, Oph N only contained a few low-mass 
high-density cores. It seems likely that the small 
star-forming molecular cores in Oph N are the sur- 
viving remnants of the GMC which produced Up- 
per Sco and are now slowly being ablated by the 
massive stars. The low star formation rate is due 
to the lack of dense cores. 
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A. Comments on individual regions 



A.l. OphN 1 

The two extinction peaks in OphN 1 lie in Lynds opacity class 6 cloud L260 (Fig. [TTj) . The YSOc and 
starless core are both in the more massive north ern extinction regio n (9 Mq above an extinction limit of 
Ay > 10, compared to a C 18 masses of 11.4M0 lNozawa et al.lll99ll ). The southern extinction peak shows 
no evidence for star formation or dense cores and little high-extinction gas (IMq above Ay — 10), indicating 
that the C 18 mass of 17.9M Q mainly lies at lower extinctions. 

The Class I source YSOcl is identified with IRAS 16442-0930/IRAS 16449-1001 and lies close to the peak 



of the northern core ( Tachihara et al 



2000b C ls O core u2). This YSOc was first identified as such from 
(| 19921) . The source is detected but unresolved in 1.3mm continuum 



the IRAS catalogue by iCarballo et al. 

with a mass < 0.05 M Q and at 850/zm by SCUBA with a peak flux less than 72 mJy/beam (L260 SMM1 



Visser et al.ll2002l see our Fig. [TU|) . c onsistent with the mm/submm emission l argely originating fr om a disk 
Ther e is an associated 12 CO outflow ( Bontemps et alJll996tlVisser et al. 2002 ) and a H2O maser ( Han et al 



1998 ). The IRAS spectrum of this source was studied in detail by iParker |l99ll) . who concluded from a 
shoulder in the short-wavelength emission at 1.5-2yLtm that it was a Class II YSO obscured by an extended 
dark cloud, but we see no evidence for this shoulder in the 2MASS fluxes. 

The submm em ission and the extinction in this region peaks towards a second core 3' to the east 
(jVisser et al.l 12002. L255 SMM2). No source is detected by Spitzer even at 70/xm (Fig. [TU]), confirming 



this as a starless core. 

A second IRAS source IRAS 16454-0955 with am biguous colours typic al of a YSO or EGO lies in the 
southeast of the IRAC region unobserved by MIPS (|Carballo et al.1 Il992h (marked in Fig. fTO)) . A search 
through the positional matches within 10" in the Spitzer catalogue found no reddened or rising spectra. 
Better waveband coverage is needed to pinpo int the source of the infrared emission detected by IRAS. 
Another ambiguous IRAS YSO or galaxy from ICarballo et al.l (|l992l) . IRAS 16439-0945, lies to the west on 
the edge of the OphN 1 A v ~ 2 region, but off the IRAC and MIPS maps. 



A.2. L234E 

South of OphN 1, L234 (Lynds opacity class 3) contains no YSOcs. The area mapp ed by Spitzer has a 
mass of only IMq above Ay > 10, though it is coincident with Tachihara et al.l (l2000b|) core t (C 18 Q mass 
I2.9M0). The SCUBA map shows weak 850/mi emission peaking at 130 m Jv/beam (jKirk et al.ll2005h . The 
faint 70^m source seen an arcminute to the north in Fig. [TOj at 16 h 48 m 06 i !6 — 10°56'06", is classified by the 
c2d pipeline as a galaxy. 

A.3. OphN 2 

This region in L204 continues the filament of bright-rimmed clouds, and contains tw o extinction peaks 
over Ay — 10 (8 and lOM respectively), several C ls O cores ( Tachihara et al. 2000bl sl-5) with masses 
from 5-33Mq, but no YSOcs. No submm maps exist of this region. It is a good candidate for future star 
formation. 



A.4. OphN 3 

Five Class II and III protostars are identified in the OphN 3 IRAC fields, four in the North and one 
in the South field. The northern sources are associated with Lynds opacity class 6 cloud L162 but are not 
coincident with the extinction peaks: they lie on the edge of the A v = 3 region heated by C Oph (see Fig.[TT|). 
One of the northern sources, Oph N YSOc3, is associate d with IRAS 16455—1405, already identified as a 
YSO or (ambi guously) an EGO by Carballo et al.l (|1992l ). This star is in fact a binary with 1" separation 
(Reipurth 44; iReipurth fc Zinnecken 19931) and possibly even a triple system ( Koreskol 12002 ). Another 
potential binary is OphN YSOc2, which lies 9" from the extended red variable T Tauri star V* V2507 Oph; 
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both have similar proper m otions consistent with the idea that these two stars form binary Reipurth 43 
(jReipurth fc Zinneckerlll993T ). 

The classical T Tauri star V* V1121 Oph (IRAS 16464-1 464) lies 10 arcminutes to the southeast of the 
other northern OphN 3 YSOcs at 16 h 49 m 15?3 - 14°22'08"6 (|Andre fc Montmerlel[l99l see our Fig. El;. It 
is detected by IRAC bands 1 and 3 but lies outside the field of IRAC 2 and 4 and the YSOc identification 
comes from the K$ — [24] selection (Table [7]) . At 24 and 70 microns it is bright but extended and confused 
with ci rrus and the catalogu e fluxe s are bandfilled (see Figs. ITOl and Hi]). V1121 Oph is classified as a Class II 
source [Andre fc Montmerlel ( 19941 ) and, like the other Class II YSOcs 4 and 5 in OphN 3, it is bright at 
70^m ( Fig. EH). 

The lVisser et all (|2002t ) SCUBA 850/im map of the L158 and L162 regions is overlai d on the 24/j,m map 



in Fig. [10] YSOc4 (Class II) is detected as a compact 850/zm source L162 SMM 1 which IVisser et all |200S 
identify with IRAS 1 6459—1411. This is a well-studied T Tauri star with a disk (jAndre fc Montmerle 
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Andrews et~aTll2009l ). YSOc5 is not detected by SCUBA at a level of 100 mJy/beam. Neither YSOc2 or 
YSOc3 lie within the SCUBA maps. The MlPS-only Class III source Ml which lies in the southeast of 
the 70/xm map in Fig. ITTlis identified as IRAS 16484-1557. The starless core in L158 detected by SCUBA 



( Visser et al . 2002) remains infrared-dark even at 70/j,m, confirming this source as starless. The nearby IRAS 
sources IRAS16445-1352, IRAS 16442-1351 and IRAS 16439-1353 appear to be cirrus. The L158 core lies 
in the most massive of the three northern extinction peaks, which have Ay > 3 masses of 5, 1 and 13Mq, 
measuring from east to west. 

The galaxy 6dFGS gJ164828.8-141437 confusingly lies among the northern OphN 3 YSOcs, and was 
initially identified by the colour-colour YSOc selection procedure but rejected when examined by eye as 
being elliptical. It appears in the 2MASS-MIPS colour-magnitude diagram (Fig. [7|) marked as a YSOc but 
clustered with the galaxies at Kg ~ 14. Its location within 20' of the other YSOcs in the north of OphN 3 
is particularly confusing (see Fig. fT0|). 

In the southern region we identify OphN YSOc5, classified Class III, with IRAS 16462—1532. This source 
is associated with the westernmost of the two extinction peaks, mass 13M Q (Ay > 10; see Fig. [TTj) . The 
southeastern extinction core is more massive, at 21Mq, and a likely site for future star formation. There are 
no submillimctrc maps of this region to determine if it contains any dense cores. The Lynds opacity 4 and 
5 dark clouds L137 and L141 lie to the north of this region, but are not associated with any YSOcs. 

A.5. OphN 4 

The OphN 4 field is associated with L1 757, which shows no extinc tion above Ay = 10, co nsistent with 
its detection in 13 CO (|Nozawa et al.lll99ll cloud C) but not C 18 (|Tachihara et all l2000bl) . There is a 



bright-rimmed globule in the 24^m map but no YSOcs. Ambiguous YSO/EGO IRAS 16355-1807 is located 
at the northern end o f the extinction filament 2 degrees east of OphN 4, at an optical extinction of Ay ~ 1 
(ICarballo et al.lll992h . 



A.6. OphN 5 

The OphN 5 field contains no mass above Ay > 10 but no YS Ocs. The dark clou d L1752 is centred 



further to the east. Op hN 5 is also associated with 13 CO cloud C (jNozawa et al.lll991l ) but no C 18 core 



(jTachihara et alJl2000b[ ) 



A.7. OphN 6 

OphN 6 is associated with the L1719 dark cloud and two high extinction regions, iTachihara et all (|2000bh 



C 18 cores a (south) and b (north), and contains two Class III YSO candidates and one flat-spectrum 
source (YSOc8). YSOc7, associated with IRAS 16186-1953, is probably a backgrou nd giant due to its 
prope r motions (see Sect. 13.51) . YSOc7 is variable in the optical/near- infrared; the IDENIS Consortium 
(2005) reports variation of more than a magnitude in / and J bands. YSOc8 is associated with IRAS 16191- 
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1936. OphN YSOc9, has a small 24/mi excess and was previously unknown. The YSOcs lie away from the 
extinction peaks, but two lie on the edges of high extinction regions: Class III YSOc7/IRAS 16186-1953 in 
the south, and the flat-spectrum YSOc8 at the edge of the northern extinction peak. 

19961) . 



The flat-spectrum source YSOc8 is associated with L1719B and has a CO outflow ([Bontemps et al 



The submillimetre e mission is still uncertain: while there is a 220 mJy dete ction at 1.3mm suggesting a disk 
mass of 0.04 Msun (jAndre &: Montmerle)ll994 : Andrews &: Wiiliamsfl2007h . SCUBA observations only give 
an 850/xm upper limit of 54 mJy/beam ( Kirk et ail 20051) . Weak submillimetre emission was detected with 
SCUBA towards L1719D at 16 h 21 m 09!2 - 20°07'10" with a flux of 62 mJy/beam. Th ese weak emission 
SCUBA jiggle maps do not reproduce well as contours and are not shown in Fig. 1101 see iKirk et al. (2005) 
for images. 

OphN 6 hosts two of the MIPS+2 MASS identified YSO cs: M3 and M4. Both of these have proper motions 
consistent with Ophiuchus distances ( Roeser et al. 2010l ). M4 lies 9" away and within the positional errors 



of ROSAT X-ray source 2RXP J162230.1-200256, supporting its identification as a YSO. 

The bright region in the southwest of OphN 6 contains four late- type B /early A stars: HIP 80019, 
HIP 80024 and HIP 80062/3. The revised Hipparcos parallaxes (jvan Leeuwenl 120071 ) and Tycho-2 proper 
motions for these stars are given in Table [IT] The weighted mean of the four parallaxes gives a distance of 
164l*g pc, larger than our assumed 130 pc (equivalent to a parallax of 7.4 milliarcsec onds) but consi stent 
with parallax distances for other stars illuminating ne bulae in the vicinity o f Sco/Oph (|Mamaiekll20081 ) and 



with distance estimates for Upper Sco (145 ± 2 pc; Ide Zeeuw et al.J 19991). The proper motions are also 



simil ar to each other and consistent with membership of Upper Sco (|Mamaiekll2008t iHoogerwerf fc Blaauw 
2000). Three of these stars are clearly associated with the cloud through nebulae (see the 3-colour images 
and 24 /j,m emission, Figs. H [3 and [JO]). HIP 80019 has no apparent nebula. HIP 80062/3, a wide B 9 binary 
with HIP 80063 itself a close binary, illuminates the optically- visible reflection nebula IC4601 (jMagakian 
2003f ). 

It is not yet clear if these stars are pre-main-sequence. None of them have the red colours characteristic 
of YSOs with disk s (Fig. [5] and CD. Nei ther HIP 8 0019 nor 80024 show emi ssion lines characteristic of 
Herbig AeBe stars ( Hernandez et al. 2005 ). However, Dahm fc Carpenter ( 2009f) found a small excess above 
10/im for HIP 80024 in Spitzer spectroscopy and characterised this source as a debris disk with the proviso 
that they 'cannot rule out that the dust is remnant primordial material'. This star lies outside the region 
mapped at 24/zm in our study so we can neither confirm or deny this disk excess. Now that their revised 
Hipparcos parallaxes and Tycho-2 proper motions are consistent with membership of Upper Sco, the two 
stars HIP 80062/63 driving IC4601 would be interesting targets for further studies. 

The bright star saturated in the IRAC bands at 16 h 20 m 20?9-20°08'05" is identified as 2MASS J16202094-2008059, 
and also saturated in 2MASS H and K, g. It is not obviously associated with any local nebulosity and, given 
its extremely red colours (V — K ~ 10. IZacharias et al.ll20051 ). is probably a background giant. 



A.8. CB68 

The CB68 globule (L146) contains a well-studied isolated Class object, from the IRAS faint source 
catalogue IRAS 16544—1604. The associated submm emission is bright and extended 40" in size with 
peak fl ux densities of 2.1 and 0.42 Jy/beam at 450 and 850 /im respectively and an estimated mass of 
0.1M© ( Huard et al. 19991 : Vallee et al. 2003 : Young et al.ll2006l ). The molecular core is elongate d northeast- 
southwest, with a perpendicular molecular outflow detected in 12 CO (|Vallee et alJ l2000L 120031 ). The dust 
emission shows strong (11%) polarisat ion yielding an estimated magnetic field strength of 120 -130 /iG in the 



radial direc tion IVallee fc Fiegel (120071) . CB68 also has a C 18 detection (jNozawa et al.lll99ll 13 CO core 34) 



though the Tachihara et al.l ( 2000bl ) C 18 Q core q2 lies 0.4 degrees further to the north. 

CB68 has more than one associated source in the Spitzer catalogue, as shown in Fig. [14] The main source 
YSOclO is strongly detected by IRAC and MIPS. A very red SED, with no 2MASS detections and rising 
towards 100^m confirms the Class status (Fig. [6]). To the southeast, 15" away, a second source is detected 
by IRAC. This lies within the point spread of MIPS emission from the central source, and is bandfilled at 
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24/im, and also shows some MUXBLEED artifacts to the north. From its SED, this source is a star (Fig. [14"]) 
and though its MIPS flux is uncertain due to the confusion with IRAS 16544—1604 there is no evidence for 
a disk excess. At 15" to the southeast, it has no associated 450/^m emission. All four IRAC bands show the 
outflow cavity extending to the northwest (Fig. [T4l . There is also nebulosity to the west and east in the 3.6 
and 4.5/im bands. 

IRAS 16534—15 57 lies fiftee n arcminutes to t he nor thwest of CB68. This source is associated with 
Lee fc Mversl l| 19991 ) core 224 or iTachihara et all (l2000bl) core q2 (see Fi g. 2 in IVallee fc Fies "' 



KM) and 

has a rising IRAS spectrum characteristic of a YSO (jCarballo et al.lll992l ). but lies beyond the range of our 
Spitzer map. 
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Fig. 1. — Regions observed by IRAC (white boxes) and MIPS (black boxes) overlaid on visual extinction from iDobashi et al.l (|2005) 
Extinction contours are at Ay = 3 and 6. OB stars associated with Upper Sco are marked by stars (|Mamaiekll2008r ) . 




Fig. 2. — Three-colour mosaics of the Oph N cloud regions imaged by Spitzer with IRAC 3.6^m (blue), IRAC 4.5/im (green) and IRAC 8.0/^m 
(red). Only the overlap regions imaged at all three wavelengths are shown. The scale bars represent 0.2 pc at the assumed distance of 
130 pc. 




Fig. 3. — Three-colour mosaics of the Oph N cloud regions imaged by Spitzer with IRAC 4.5/^m (blue), IRAC 8.0/im (green) and MIPS 24/^m 
(red). Only the overlap regions imaged at all three wavelengths are shown. The scale bars represent 0.2 pc at the assumed distance of 
130 pc. 
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Fig. 4. — Differential source counts compared to the Wainscoat model for the galactic stellar background. 
The thick grey line shows the differential source counts per square degree in each band for all catalogued 
sources. The thin black line shows the differential source count for only the sources classified as stars. The 
stellar count can be compared to the dashed line which shows the Wainscoat model for the stellar density 
in the Galaxy. 
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Fig. 5. — Co lour-magnitude diag rams for all Oph N regi ons combined, showing the YSOc identification 
criteria from Harvey et al.l (|2007f ) and lEvans et al.l ( 2007b . Symbols are pointlike/extended YSOcS (red 
circles/crosses) and pointlike/extended galaxies (green/purple) as shown in the key. The black lines show 
the fuzzy colour-magnitude cuts that define the YSO candidate criterion. Objects fainter than the hard 
limits shown as red dot-dashed lines are excluded from the YSOc category. The area covered is the 1.68deg 2 
covered by both IRAC and MIPS. 
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Fig. 6.— Spectral energy distributions for the ten YSOcs listed in Table El Measured NOMAD, 2MASS, 
IRAC and MIPS fluxes are shown as open circles. For the evolved sources (flat-spectrum, Class II and Class 
III), photospheric models are plotted for a K7 star except for Sources 3 and 4, for which an AO photosphere 
was used. The filled circles (green in the online version) then show the fluxes corrected for reddening. For 
the Class I sources 1 and 8, no photospheric model is plotted. The dashed line sh ows the median SEP o f 
T Tauri stars in Taurus (with error bars showing the quartiles of the distribution) ( Hartmann et al. 2005f) . 
Uncertainties are smaller than the symbol size (see Table 
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Fig. 7. — Ks vs. Ks — [24] colour-magnitude diagrams for all Oph N regions combined. Symbols reflect the 
c2d pipeline classification: stars (grey dots); star+dust (red plus symbols in the online version); galaxies (x, 
green in the online version); YSOc (black circles, filled red in the online version); and other (black dots). 
Class III, Class II, fiat-spectrum and Class I source models and the horizontal and diagonal lines mark cuts 
for faint Ks and 24/mi detections. 
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Fig. 8. — Spectral energy distributions for sources detected via K-24/jxn colours (see Table [7|)- Fluxes from 
2MASS and Spitzer are plotted as open circles and those from WISE as filled circles (red in the online 
version). 
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Fig. 9. — Colour-magnitude diagrams for the IRAC-detected YSOc (open circles) comparing the AGB 
candidates (filled c ircles) with the spec troscopically identified AGB stars in Serpens (crosses, Harvey 2012 
priv. comm., from loiiveira et al.ll2009h . The color-magnitude combinations are the same as Fig. [5j though 
unlike that figure the galaxy interloper 6dFGS gJ164828.8 -141437 is not shown. The YSOc occupying the 
same region of the graph as the AGB stars, but with the lowest 24^m fluxes and bluest [4.5] — [8.0] colours, 
is YSOc 9. 
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Fig. 10. — SCUBA 850/^m contours (green in online version, from 100 mJy/beam in steps of 100 mJy/beam) 
overlaid on MIPS 70/im maps (grey/colourscale, logarithmic from —10 to 5 MJy/sterad). Top to bottom, 
left to right: (a) OphN 1, (6) OphN 3b, (c) OphN 6, (d) OphN 3a, (e) L234E and (f) CB68. The regions 
mapped by IRAC (band 1) are marked in black. Circles and triangles represent IRAC-detected and MIPS- 
only YSO candidates from Tables [5] and [7] respectively. YSOcs identified using IRAC and MIPS are marked 
with circles, and those identified using MIPS and 2MASS with triangles (in the online version, coloured 
according to YSO class: red=Class I, yellow=flat-spectrum; pink=Class II; blue=Class III). The T Tauri 
star V1121 Oph, the confusing galaxy gJ164828.8 -141437, and the possible YSO IRAS 16454-0955 are also 
marked. The three regions which are starless and unobserved by SCUBA are not shown (OphN 2,4,5). 
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Fig. 11. — MIPS 24/im maps (colourscale) overlaid with Spitzer-derived extinction contours at Ay — 
[6,10,20] (black) for the areas observed by IRAC. Top left to bottom right: (a-f) OphN 1,2,3,4 5 and 
6, (g) L234E and (h) CB68. For compari son, the Ay = 3 level from optical starcounts ( Dobashi et al.ll2~005l ) 
and Ay = 4 from near-infrared colours ( Rowles fc Froebrichl 200!); are plotted with dashed contours (blue 
and red respectively in the online version). Sources are marked as in Fig. [TU1 
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Fig. 12. — Spitzer 160/im emission and IRAS HIRES-processed 100/im emission. Left to right, top to bottom: 
(a) IRAS 100/im image of OphN 1,2, and 3, with arrows indicating the approximate direction of C Oph; (6) 
Spitzer 160/im image of OphN 1,2 and 3; (c) IRAS 100/im image of OphN 4,5 and 6; (d) Spitzer 160/im 
image of OphN 4,5 and 6. Greyscales are a square-root stretch from -10 to 50 MJy/steradian (IRAS 100/im) 
and 10 to 150 MJy/steradian (Spitzer 160/im). Sources are marked as in Fig. [TUT 
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Fig. 13. — Excess 160/im over 100/im emission (greyscale, linear from -50 to 50 mJy/steradian) overlaid with 
Spitzer Ay — 6 contour (green in the online version, shown in the IRAC regions only). Left to right, (a-f), 
regions OphN 1, 2, 3, 4, 5 and 6. CB68 and L234E have no 160/im map and are not shown. Sources are 
marked as in Fig. [TO] 
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Fig. 14. — Top: CB68 imaged by IRAC 2 (4.5/im) showing the reflection nebula extending towards the 
northwest in the direction of the outflow. The black contours are MIPS 70/im[50,100,500] MJy/sterad. 
Markers show the YSOc (circle, red in the online version) and southern star (cross, green in the online 
version) from the Spitzer catalogue. Bottom: CB68 SEDs for the YSO (circles) and the confusing star to 
the southeast (crosses). 
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Table 1 

AORs in the Spitzer Gould Belt Oph N (Scorpius) Catalogue 



Instrument RA(J2000) Doc(J2000) PID Duration(mins) Observation date/time AOR key 



OphN 6 



IRAC 


16:21:24.19 


-19:56:01.80 


30574 


67.19 


2007-09-09 


03:51:06.6 


19992320 


IRAC 


16:21:24.48 


-19:56:05.80 


30574 


67.18 


2007-09-09 


07:49:28.5 


19991808 


MIPS 


16:21:29.60 


-19:57:45.00 


30574 


57.96 


2007-04-08 


18:13:07.4 


20001280 


MIPS 


16:21:29.60 


-19:57:45.00 


30574 


57.96 


2007-04-09 


00:39:24.0 


20000768 


OphN 5 
















IRAC 


16:27:45.78 


-19:24:43.00 


30574 


25.38 


2006-09-21 


13:19:07.2 


19960064 


IRAC 


16:27:45.78 


-19:24:43.00 


30574 


25.38 


2006-09-21 


17:17:44.2 


19959552 


MIPS 


16:27:40.00 


-19:25:22.00 


30574 


30.83 


2007-04-08 


10:55:48.0 


19968512 


MIPS 


16:27:40.00 


-19:25:22.00 


30574 


30.83 


2007-04-08 


15:39:59.8 


19968000 


OphN 4 
















IRAC 


16:31:53.68 


-19:35:03.00 


30574 


54.31 


2006-09-21 


12:12:25.1 


19990784 


IRAC 


16:31:53.97 


-19:35:07.00 


30574 


54.32 


2006-09-21 


16:27:39.8 


19990272 


MIPS 


16:32:21.50 


-19:43:36.10 


30574 


44.37 


2007-04-06 


23:54:22.0 


19999744 


MIPS 


16:32:21.50 


-19:43:36.10 


30574 


44.37 


2007-04-07 


05:44:06.0 


19998976 


OphN 3 (including c2d 


L158) 












IRAC 


16:50: 6.45 


-15:29:48.50 


30574 


68.54 


2006-09-26 


15:36:32.5 


19959296 


IRAC 


16:50: 6.83 


-15:29:47.80 


30574 


68.54 


2006-09-26 


22:11:19.0 


19959040 


MIPS 


16:49:60.00 


-15:30: 9.70 


30574 


78.19 


2007-04-07 


20:15:50.3 


19966720 


MIPS 


16:49:60.00 


-15:30: 9.70 


30574 


78.19 


2007-04-08 


04:12:25.9 


19966208 


IRAC 


16:49: 4.02 


-14:11:09.10 


30574 


48.44 


2006-09-26 


20:52:06.4 


19990016 


IRAC 


16:49: 4.40 


-14:11:08.40 


30574 


48.44 


2006-09-27 


00:51:11.2 


19989504 


MIPS 


16:48:53.60 


-14: 5:45.00 


30574 


78.19 


2007-04-08 


19:09:17.5 


19996160 


MIPS 


16:48:53.60 


-14: 5:45.00 


30574 


78.19 


2007-04-09 


01:35:34.1 


19995904 


IRAC 


16:47:08.00 


-13:58: 5.00 


139 


13.81 


2004-08-16 


18:55:38.5 


5126400 


IRAC 


16:47:08.00 


-13:57:55.00 


139 


13.81 


2004-08-17 


07:35:56.5 


5126912 


MIPS 


16:47:15.60 


-13:59:28.80 


139 


20.46 


2005-03-10 


02:36:55.6 


9415680 


MIPS 


16:47:16.40 


-13:56:41.20 


139 


20.46 


2005-03-10 


09:47:39.7 


9427456 


OphN 2 (including c2d 


L204C) 












IRAC 


16:47:32.48 


-12:06:21.40 


30574 


32.86 


2006-09-23 


01:45:10.0 


19958784 


IRAC 


16:47:32.85 


-12:06:20.60 


30574 


32.87 


2006-09-23 


06:07:36.2 


19958528 


MIPS 


16:47:33.80 


-12:09:06.40 


30574 


30.82 


2007-04-08 


10:26:38.9 


19964160 


MIPS 


16:47:33.80 


-12:09:06.40 


30574 


30.82 


2007-04-08 


15:10:50.8 


19963904 


IRAC 


16:47:39.00 


-12:21:09.00 


139 


8.49 


2004-08-17 


07:47:01.5 


5127424 


IRAC 


16:47:39.00 


-12:20:59.00 


139 


8.49 


2004-08-18 


05:38:13.0 


5127936 


MIPS 


16:47:42.30 


-12:21:11.70 


139 


20.39 


2005-03-10 


02:18:26.2 


9408512 


MIPS 


16:47:35.70 


-12:21:06.30 


139 


20.39 


2005-03-10 


09:29:20.8 


9422080 


OphN 1 
















IRAC 


16:47:24.26 


-9:46:39.40 


30574 


43.07 


2006-09-23 


18:42:25.8 


19988224 


IRAC 


16:47:24.63 


-9:46:38.60 


30574 


43.07 


2006-09-24 


01:31:23.5 


19987712 


MIPS 


16:47:24.01 


-9:52:44.60 


30574 


30.82 


2007-04-08 


17:43:41.6 


19994368 


MIPS 


16:47:24.01 


-9:52:44.60 


30574 


30.82 


2007-04-08 


23:11:15.0 


19994112 


L234E 
















IRAC 


16:48:09.00 


-10:55:56.00 


139 


12.05 


2004-08-17 


07:53:06.1 


5128448 


IRAC 


16:48:09.00 


-10:55:46.00 


139 


12.04 


2004-08-18 


05:44:20.5 


5128960 


IRAC 


16:48:12.20 


-10:55:59.60 


139 


20.39 


2005-03-10 


02:00:09.4 


9411072 


IRAC 


16:48:05.80 


-10:55:52.40 


139 


20.39 


2005-03-10 


09:10:57.7 


9431808 


L146 / CB68 














IRAC 


16:57:20.50 


-16:09:02.00 


139 


8.83 


2004-08-17 


08:23:08.3 


5130752 


IRAC 


16:57:20.50 


-16:08:52.00 


139 


8.83 


2004-08-18 


07:55:40.2 


5131264 


MIPS 


16:57:20.10 


-16:10:25.80 


139 


18.87 


2005-04-04 


05:55:18.0 


9425408 


MIPS 


16:57:20.90 


-16:07:38.20 


139 


18.88 


2005-04-05 


22:55:58.5 


9439232 



Note. — The original AOR names of regions OphN 1—6 were Sco 1—6. 
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Table 2 

Associations with Lynds Dark Nebulae and molecular cores. 



Region CO a LDN (opacity) b 

OphN 1 U, ul,u2 L260 (6), L255 (6) 

L234E T, t L234 (3) 

OphN 2 S, sl-5 L204 (6) 

OphN 3a R, rl,r2 L152 (3), L158 (6), 

L162(6), L163(2), L141 

(5),L137 (4) 

OphN 3b P, pl,p2 
CB68 Q, ql L146 (5) 

OphN 4 C - L1757 (6) 

OphN 5 C,- L1752 (4) c 

OphN 6 A, B, a, b L1719 (5) c 



a Co res identified from 13 CO/C 18 (jTachihara et al 
2000bh 



3 Lynds Dark Nebula (|Lvndslll962l ) 



c Dark Nebulae L1719 and L1752 lie in the constel- 
lation of Scorpius, whereas all the others listed here 
are within Ophiuchus. 
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Table 3 

Detection statistics with IRAC and 2MASS. 



Detection type Detections 



At least one IRAC band: 106200 

All 4 bands : 5378 

Only 3 bands : 7902 

Only 2 bands : 50632 

Only 1 bands : 42288 

MIPS 1 845 

MIPS 1 and IRAC 762 

MIPS and all 4 IRAC bands 449 

2MASS k 10919 

2MASS k and h 10822 
2MASS k and at least one IRAC band 10917 

2MASS k and all 4 IRAC bands 4723 

2MASS k but not IRAC 2 

IRAC band 1 but not 2MASS k 89038 



MIPS, all 4 IRAC bands, and 2MASS k 196 



Note. — Detections are counted above a signal-to-noisc 
ratio of 3 and in the overlap area of all 4 IRAC bands 
only. 
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Table 4 

Detection of sources above various S/N thresholds in individual IRAC bands. 



S/N 


3.6^m 


4.5/xm 


5.8^m 


8.0/im 


S/N of 3 


136853 


94166 


18320 


821 


S/N of 5 


101913 


64178 


9618 


5465 


S/N of 10 


50565 


28830 


4742 


2696 


S/N of 15 


23182 


13459 


2446 


1460 
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Table 5 

yso candidates with irac detections 



YSOc 




Field 


Class 


a 


ID (offset) 


S u. t) ii "i in core 


I^C f C V C 11 c c s 


1 


J164658. 3-093519 


OphNl 




0.66 


IRAS 16442-0930 (0.9') 


L260 SMM1 


1,2,3 


2 


J164817. 6-141109 


OphN3 




-0.70 


V2507 Oph(9.0"), Oph4 (9.0") b 




4,5 


3 


J164821. 9-141043 


OphN3 




-1.03 


IRAS 16455-1405 (0.27"), PDS 91 




1,6 


4 


J164845. 6-141636 


OphN3 




-0.97 


IRAS 16459-1411, V2508 Oph, Oph6 (0.16") 


L162 SMM1 


5,7,8 


5 


J164900. 8-141711 


OphN3 




-1.28 


Oph 5 (0.12"), HBC 654 




5 


6° 


J164905. 6-153713 


OphN3 


III 


-1.87 


IRAS 16462-1532 (7.45") 






7° 


J162137. 7-200037 


OphN6 


III 


-2.09 








8 


J162204. 3-194327 


OphN6 


F 


0.02 


IRAS 16191-1936 (1.04") 


L1719B 


1,2,9 


9 


J162209. 6-195301 


OphN6 


III 


-2.11 








10 


J165719. 6-160923 


CB68 




2.40 


IRAS 16544-1604 (12.05") 


CB68,CB68SMM1 


10,1,11 



a Infrared spectral index a from a fit to fluxes from K$ band to MIPS 24/j.m (Evans et al.l f2007). 
b The Spitzer detection and V2507 Oph arc separate stars, probably in a binary. Sec Appendix. IA.41 
c Probably an AGB star. See Sect. [33] 



Note.— References : (1) Carballo c <~aTl 1119921); (2 ) Bontcmps ct a T] fll996Tk (3) IVisser et al.l IpOOll) ; (4) IDucourant et all 120051) ; (5) 
Andre fc Montmerlej i 199411 i ; (6) IReipurth fc Zinneckeri U993I) : ( 71 1 Andrews et al] 120091) ; (81 IVisser et al.l t5002T) ; f9) IKirk et al.l 12003) ; (10) 
Clemens fe Barvainisl 1 119881) ; qi lHuard et al.l | |1999|K 



Table 6 

Fluxes for YSO candidates with IRAC detections. 



YSOc 




2MASS 






Snifacr TR AC! 






Spitzer MIPS 
F70 
mjy 




J 

mjy 


H 

m.Ty 


K s 
mjy 


F3.6 
m.Ty 


F4.5 
mjy 


F5.8 
mjy 


F8.0 
mjy 


F24 
mjy 


F160 
mjy 


1 


1.8 ± 0.1 


8.3 ± 0.3 


26.8 ± 0.8 


81. 8± 4.3 


150. 0± 7.8 


231± 11 


189± 12 


2430± 248 




3700± 730 


2 


50.9 ± 1.7 


90.4 ± 4.4 


105.0 ± 2.6 


71. 4± 3.5 


64. 6± 3.1 


68± 3 


90± 4 


168± 16 


87± 12 


< 2200 


3 


123.0 ± 3.2 


298.0 ± 7.7 


477.0 ± 7.9 


641.0±31.9 


559.0±30.5 


539± 26 


471± 23 


639± 59 


476± 49 


< 2200 


4 


519.0 ± 11.5 


957.0 ± 33.5 


1200.0 ± 29.9 


1010.0±51.5 


889.0±46.0 


816± 40 


925± 44 


1310± 124 


873± 92 


< 2200 


5 


371.0 ± 8.2 


571.0 ± 21.0 


575.0 ± 11.1 


390.0±19.9 


315.0±19.1 


274± 13 


270± 13 


303± 28 


92± 14 


< 2200 


6 


140.0 ± 3.2 


663.0±25.6 


1120.0±21.8 


1000.0±51.6 


778.0±43.4 


742± 36 


517± 25 


146± 14 


< 11 


< 2200 


7 


437.0± 8.4 


1150.0±38.2 


1510.0±25.0 


1100.0±57.9 


809.0±40.4 


694± 34 


472± 22 


114± 11 


< 11 


< 2200 


8 


27.8± 0.7 


136.0± 2.9 


347.0± 8.3 


759.0±38.3 


1090.0±56.3 


1440± 69 


2010± 110 


3730± 390 


5320± 573 


8700± 730 


9 


144.0± 3.8 


208.0± 4.8 


184. 0± 4.1 


109. 0± 6.6 


71. 5± 4.1 


50± 3 


33± 2 


16± 2 


< 11 


< 2200 


10 











0.4± 0.0 


1.4± 0.1 


1.9±0.1 


1.1±0.1 


659±61 


13100±1210 





a Sourcc lies outside the 70//m map. 



Table 7 

yso candidates identified from mips 24 micron and 2mass in regions with no irac coverage. 



YSOc 


SSTgbs 


Field 


Class 


ID(offsot) 


References 


Ml a 


J165121. 4-160259 


OphN3 


III 


IRAS 16484-1557 (7"), ASAS 165122-1602.9 (5.9") 


1 


M2 a 


J163142. 5-203814 


OphN4 


II 






M3 


J162214. 8-204540 


OphN6 


II 






M4 


J162229. 8-200248 


OphN6 


II 


2RXP J162230. 1-200256 (9.1") 




M5 b 


J164915. 3-142209 


OphN3 


II 


V1121 Oph(0.08") 


2.3 



a Probably an AGB star. Sec Scct.l3~5l 

b V1121 Oph lies outside the IRAC mapping region at 4.5 and 8.0^tm so is not classified as a YSOc by the c2d pipeline. 



Note.— References: (l) IPoimanski et al.l (120051'! : (2) IAndre fe Montmerle] 1119941) : fS lCarballo et alj J1992TI . 



Table 8 

Fluxes for YSO candidates identified from MIPS 24 micron and 2MASS. 



YSOc 




9MASS 








WTRF, 






Spitzer MIPS 
F70 
mJy 




J 

mjy 


H 

mjy 


Ks 
mjy 


F3.4 
mjy 


F4.6 
mjy 


F12 
mjy 


F22 
mjy 


F24 
mjy 


F160 
mjy 


Ml 


969.0 ± 18.7 


1550.0 ± 28.5 


1750.0 ± 33.9 


1029.7 ± 38.9 


1086.9 ± 24.0 


406.91 ± 7.12 


185.34 ± 4.27 


178 ± 17 


< 11 


< 2200 


M2 


93.1 ± 1.8 


203.0 ± 3.9 


235.0 ± 4.6 


202.6 ± 2.4 


154.7 ± 2.4 


158.31 ± 2.62 


93.06 ± 2.40 


77 ± 7 


40 ± 7 


< 2200 


M3 


21.8 ± 0.5 


24.0 ± 0.6 


21.6 ± 0.4 


11.8 ± 0.3 


8.5 ± 0.2 


4.79 ± 0.21 


6.92 ± 1.21 


7± 1 


< 11 


< 2200 


M4 


115.0 ± 2.3 


215.0 ± 4.2 


243.0 ± 4.3 


222.2 ± 4.9 


194.2 ± 3.6 


115.75 ± 1.92 


128.10 ± 2.95 


112 ± 10 


120 ± 16 


< 2200 


M5 


762.0 ± 27.4 


1070.0 ± 23.6 


1100.0 ± 26.2 


746.6 ± 22.7 


744.4 ± 15.1 


3170.32 ± 73.00 


5915.10 ± 49.03 


3850 ± 221 a 


2400 ± 79 a 


2050 ± 730 



a 24 and 70//m fluxes from 100" radius aperture photometry. The listed c2d pipeline point-source 24/im flux was 410 ± 200 mjy. 



Table 9 

Core masses above extinction thresholds. 



Region 


M{A V > 6) a 


M(A V > 10) 


M(A V > 20) 




M G 


M Q 


M Q 


OphN 1 


30 


10 


1 


OphN 2 


51 


17 


2 


OphN 3 


149 


52 


7 


OphN 4 


9 








OphN 5 


G 








OphN 6 


40 


5 





CB68 


1 








L234E 


8 


1 







294 


85 


10 



a The Ay — 6 masses are lower limits as these regions 
extend beyond the limits of the map. See Fig. [11] 
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Table 10 

100/160 micron spectral slope and average dust temperatures 



Region 


slope a 


T d /K b /0 = 2.0 


OphN 1 


0.35 


15.6 


OphN 2 


0.39 


16.1 


OphN 3S 


0.37 


15.9 


OphN 3N 


0.40 


16.3 


OphN 4 


0.39 


16.1 


OphN 5 


0.45 


17.0 


OphN 6 


0.44 


16.8 



Uncertainty in slope is ±0.002. Esti- 
mated systematic uncertainty of 10 per- 
cent leads to ±0.4 K absolute tempera- 
ture uncertainty. 

b Color-corrected dust temperature for 
= 2.0. Changing from = 
2.0 to = 1.7 (appr opriate for 



Ossenkopf &: Henningl 119941 model OH5 



dust opacity) leads to dust temperatures 
that are 0.8 K systematically higher. 
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Table 11 

Parallaxes and proper motions for early- type stars in OphN 6. 



HIP 


Spectral type 


Parallax 


mu a 


US 






mas 


arcsec/yr 


arcsec / yr 


80019 


B9 a 


5.17 ±1.10 


-12.01 ±1.38 


-25.54 ±1.01 


80024 


B8 b 


6.12 ±0.78 


-11.89 ±1.01 


-24.25 ±0.71 


80062 


B9.5 Va c 


9.82 ± 2.44 


-17.31 ±2.88 


-19.40 ±2.08 


80063 


A0 d 


7.45 ± 3.42 


-17.55 ±4.17 


-13.24 ±2.86 



No te. — Parallaxes and proper motions are taken from Ivan Leeuwen 
(|2007l) . 

^Hernandez et al.l (|2005l ) 

b Strictly, kB8HeA0hB8IImA2IbSi(Cr:) (jGarrison fc Gravlll994l ). 
C HIP 80062/63 form a visual double (|Grav et al.ll2006h 
d kA0hA0VmA2 (|Abt et al.lll979h 
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